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ABSTRACT 

Three periodically variable stars have recently been discovered (V407 Vul, P = 
9.5 min; ES Cet, P = 10.3 min; RX J0806. 3+1527, P = 5.3 min) with properties 
that suggest that their photometric periods are also their orbital periods, making 
them the most compact binary stars known. If true, this might indicate that close, 
detached, double white dwarfs are able to survive the onset of mass transfer caused 
by gravitational wave radiation and emerge as the semi-detached, hydrogen-deficient 
stars known as the AM CVn stars. The accreting white dwarfs in such systems are 
large compared to the orbital separations. This has two effects: first it makes it likely 
that the mass transfer stream can hit the accretor directly, and second it causes a loss 
of angular momentum from the orbit which can destabilise the mass transfer unless 
the angular momentum lost to the accretor can be transferred back to the orbit. The 
effect of the destabilisation is to reduce the number of systems which survive mass 
transfer by as much as one hundred-fold. In this paper we analyse this destabilisation 
and the stabilising effect of a dissipative torque between the accretor and the binary 
orbit. We obtain analytic criteria for the stability of both disc- fed and direct impact 
accretion, and carry out numerical integrations to assess the importance of secondary 
effects, the chief one being that otherwise stable systems can exceed the Eddington 
accretion rate. We show that to have any effect upon survival rates, the synchronising 
torque must act on a timescale of order 1000 years or less. If synchronisation torques 
are this strong, then they will play a significant role in the spin rates of white dwarfs 
in cataclysmic variable stars as well. 

Key words: binaries: close — accretion, accretion discs — gravitational waves — 
white dwarfs — novae, cataclysmic variables 



1 INTRODUCTION 

In recent years there have been discoveries of close detached 
double white dwarfs at a rate that suggests that there may 
be of order 200 million such systems in our Galaxy, making 
them the largest population of close binary stars (e.g. Marsh, 
Dhillon & Duck 1995; Napiwotzki et al. 2003). These sys- 
tems were predicted theoretically long ago (Webbink 1979; 
Webbink 1984; Iben & Tutukov 1984), and were proposed 
as potential Type la supernova progenitors, the explosion 
being triggered by the action of gravitational wave losses 
which cause them to merge. Only those systems with to- 
tal mass in excess of 1.4 Mq have the potential to become 
Type la supernovae. These represent only a few percent of 
merging systems (Iben, Tutukov & Yungelson 1997), which 
raises the question as to the fate of the less massive systems. 
Possible outcomes which have been suggested are sdB stars 
and R CrB stars (Webbink 1984; Iben 1990; Saio & Jeffery 



2002). Another possibility is that the stars survive as binary 
systems to become the semi-detached accreting binary stars 
known as the AM CVn stars (Nather, Robinson & Stover 
1981; Tutukov & Yungelson 1996; Nelemans et al. 2001). 
These remarkable systems have extremely short orbital pe- 
riods (17 to 65 minutes) made possible because of the high 
densities of their degenerate donor stars. Moreover they fea- 
ture accretion discs composed largely of helium which offer 
the opportunity to see the effect of abundance upon disc 
physics. 

In their study of the evolution of AM CVn systems, 
Nelemans et al. (2001) considered the evolution of detached 
double white dwarfs into the semi-detached phase, which re- 
quires the systems to pass through an ultra-compact phase 
when orbital periods as low as two or three minutes are pos- 
sible. The stability of the mass transfer during this stage is 
crucial as to whether the systems survive intact or merge. 
An additional complication is that even in stable systems 
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the equilibrium mass transfer rate can be (highly) super- 
Eddington, resulting in mass and angular momentum loss 
and probably merging of the system (Han & Webbink 1999; 
Nelemans et al. 2001). Mass transfer between two white 
dwarfs has generally been taken to be stable if the mass 
ratio q = M2/M1, where Mi is the mass of the accretor, 
is smaller than about 2/3. However, Nelemans et al. iden- 
tified the interesting possibility of the mass transfer stream 
hitting the accreting white dwarf directly during this phase. 
They realised that mass transfer without an accretion disc 
may destabilise the mass transfer process for lower mass ra- 
tios still because the normal transfer of angular momentum 
from disc to orbit would be absent. 

Assuming that unstable systems do not survive, Nele- 
mans et al. (2001) found that the destabilisation by direct 
impact accretion could have a dramatic effect upon the num- 
ber of systems able to survive the onset of mass transfer. If 
close white dwarfs typically do not survive this phase, then 
it may mean that either evolution from helium stars (Iben 
& Tutukov 1991; Nelemans et al. 2001) or from immediate 
post main-sequence donors (Podsiadlowski, Han & Rappa- 
port 2003) are more important channels for the formation 
of AM CVn stars than the double white dwarf route. 

There is now some observational evidence, that the 
very short orbital periods suggested by the white dwarf 
merger model are possible. The three stars V407 Vul (Crop- 
per et al. 1998), ES Cet (Warner & Woudt 2002), and 
RX J0806.3+1527 (Ramsay, Hakala & Cropper 2002; Israel 
et al. 2002) show periodic variations on periods of 9.5, 10.3 
and 5.3 min respectively, and in each case no other period 
is seen. In addition, there is evidence that each of these 
stars appears hydrogen deficient (strongest for ES Cet), as 
expected if these periods are truly their orbital periods. Fi- 
nally, V407 Vul and RX J0806. 3+1527 show several obser- 
vational characteristics compatible with their still being in 
a phase of direct impact accretion (Marsh & Steeghs 2002; 
Ramsay et al. 2002; Israel et al. 2002), or alternatively that 
they are still detached but approaching mass transfer (Wu 
et al. 2002; Strohmayer 2002; Strohmayer 2003; Hakala et al. 
2003). What is still lacking is spectroscopic proof of the 
ultra-short periods of these stars, nevertheless it is likely 
that to understand these systems we will need to understand 
mass transfer between white dwarfs. 

In this paper we study the stability of mass transfer 
in double white dwarf binary stars. In particular we study 
whether, and under what circumstances, dissipative cou- 
pling of the accretor to the orbit, through either a magnetic 
field or tidal forces, can stabilise the mass transfer. We find 
that the finite size of the accretor not only leads to direct 
impact accretion, but also destabilises the accretion whether 
it occurs by direct impact or through a disc. The destabil- 
isation has a significant effect upon the survival of double 
white dwarfs as binary stars. 

We start in section |2] by setting up the equations de- 
scribing the onset of mass transfer. Then, in sections[3]and2] 
we present analytical and numerical solutions of these equa- 
tions. In section|S]we show the possible impact of our find- 
ings for the population of AM CVn systems while in sec- 
tion|S|we discuss the uncertainties and open questions. 



2 THE EQUATIONS GOVERNING THE 

EVOLUTION OF MASS TRANSFER AND 
WHITE DWARF SPIN 

To consider the stability of mass transfer, we need a model 
that includes a description for the variation of mass trans- 
fer rate with the degree of over-filling of the Roche lobe by 
the donor star. We discuss this in section We also need 
to consider the question of feedback of angular momentum 
from the accreting white dwarf to the orbit (similar to Pried- 
horsky & Verbunt (1988) who considered feedback from a 
disc). Before mass transfer starts, the two white dwarfs will 
typically have been orbiting one another in ever decreasing 
circles for many millions of years. Tidal or magnetic cou- 
pling will be acting to synchronise their spins with their or- 
bit. We assume that the donor is always synchronised. Once 
mass transfer starts, the accretor will start to be spun up 
by the addition of high angular momentum material. What 
matters for stability is whether it is able to couple strongly 
enough to the orbit for the added angular momentum to be 
fed back into the orbit. This will be the new element that 
we add to the formalism of Priedhorsky & Verbunt (1988), 
along with a proper account of the angular momentum of 
material in the inner disc in the disc-fed case. 

To understand how the mass transfer rate changes with 
time, we first need to know how the orbital separation 
changes. This evolves because of orbital angular momen- 
tum loss, and thus we start by considering the evolution of 
angular momentum. 



2.1 Angular momentum loss 

The orbital angular momentum changes because of the ac- 
tion of gravitational radiation, the loss of mass during mass 
transfer and the coupling between the accretor's spin and 
the orbit. These can be written as: 

Jorb = Jgr + VGMiR h M 2 + kMlR * u. (1) 

TS 

The first term on the right-hand side represents the change 
from gravitational wave radiation. It is given by 

f 32 G 3 M1M2M 

JGR = -T— a* Jorb ' (2) 

where Mi and M2 are the masses of the two stars, a is the 
orbital separation and M — M1 + M2 is the total mass of the 
two stars (Landau & Lifshitz 1975). The second term on the 
right of Eq.Qrepresents the loss of angular momentum from 
the orbit upon mass loss at a rate of M2. We adopt Verbunt 
& Rappaport's (1988) description of the angular momen- 
tum of the transferred matter in terms of the radius Rh of 
the orbit around the accretor which has the same specific 
angular momentum as the transferred mass. We use Eq. 13 
of Verbunt & Rappaport (1988) to calculate Rh, accounting 
for their inverted definition of mass ratio compared to ours. 

The third term of Eq. is a torque from dissipa- 
tive coupling, tidal or magnetic, which we parameterise in 
terms of the synchronisation timescale ts of the accretor, 
with the torque linearly proportional to the difference be- 
tween the accretor's spin and the orbital angular frequencies, 
uj — Q, B — 0, o . The term kM\R\ is the moment of inertia of 
the accretor, where k « 0.2 (Motz 1952); a more accurate 
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approximation will be derived later. We ignore the spin an- 
gular momentum of the synchronised donor. 
The orbital angular momentum is given by 



Jen 



/Ga\ 1/2 
\MJ 



Mi Ms, 



from which it follows that 

Jorb / 1 \ M2 la 

Jorb" 1 q) M 2 + 2a' 



(3) 



(1) 



for conservative mass transfer (M = 0) which we assume to 
be the case throughout this paper. 

Combining Eqs [T] H and 01 one obtains 



J orb 



Jgr . rr, ; : M 2 . kMiR\ 
-j— + \/(l + q)r h — + — u), 

Jorb M 2 TS J orb 



(5) 



where rn = Rh/a, and then Eqs [I] and [S] can be combined 
to obtain an expression for the rate of change of the orbital 
separation: 



a J G r , kM 1 Rf 
0~ = 1 1 T — w 

LCI <-*orb ' S '-'orb 



-(l-c-V(IT^K)^. (6) 

This equation shows that the orbital separation decreases 
due to gravitational wave losses, increases due to the dis- 
sipative torque if the accretor is spinning faster than the 
binary, and may decrease or increase as a result of mass 
transfer, depending upon the mass ratio. The last term is 
the reason why mass transfer can be unstable because it 
can lead to a runaway where mass transfer causes the sep- 
aration to decrease which in turn causes the mass transfer 
rate to increase. To see this in detail, we now specify how 
the mass transfer rate depends upon the degree to which 
the Roche lobe is over-filled, and how the degree of overfill 
evolves with time. 



2.2 The mass transfer rate and the evolution of 
the over-fill factor 

We define the over-fill of the Roche lobe to be 



A = R 2 -R L , 



(7) 



where R 2 is the radius of the donor and Rl the radius of 
its Roche lobe. The mass loss rate from the donor mono- 
tonically increases with A. How it does so has been inves- 
tigated by many authors (Paczyriski & Sienkiewicz 1972; 
Webbink 1977; Savonije 1978; Meyer & Meyer-Hofmeister 
1983; Ritter 1988; Priedhorsky & Verbunt 1988). These in- 
vestigations divide into two classes which we will term "adi- 
abatic" and "isothermal". At low mass transfer rates, the 
inner Lagrangian point lies close to the photosphere of the 
donor, and we may adopt the formalism of Ritter (1988) to 
give 

»V A/H 



M 2 = 



(8) 



where Mo depends upon the density in the photosphere and 
H, the scale-height, is given by H — ksT 2 / 1 ^g 2 where fcg is 
Boltzmann's constant, and T 2 and g 2 are the surface temper- 
ature and gravity of the donor. We term this the "isother- 
mal" model. It is appropriate to cataclysmic variable stars 
and to AM CVn stars at long orbital periods. At high mass 
transfer rates, the inner Lagrangian point will move below 



the photosphere, and the exponential relation of the isother- 
mal case will break down. The correct way to treat this is 
to compute the stellar structure (Savonije 1978). However 
with all calculations overshadowed by the uncertainty in the 
synchronisation torque, we prefer to move to the opposite 
extreme of the isothermal case, which is the adiabatic re- 
sponse, which for a white dwarf donor gives a mass transfer 
rate of 



M 2 



-f(Mi,M 2 ,a,R 2 )A 3 



(9) 



for A > 0, and zero for A < (Webbink 1984). Combin- 
ing results from Webbink (1984), Chandrasekhar (1967) and 
Webbink (1977), the function / is given by 

f(Mi,M 2 ,a,R 2 ) = —(-7^1 (10) 



3/2 



l f 3fj,M 2 y- 1/2 

X P{^R 2 ) Ma 2 -1)] , 

where m e is the mass of an electron, m n is the mass of a 
nucleon, fi' e is the mean number of nucleons per free electron 
in the outer layers of the donor (which we will assume to be 
2), P is the orbital period, vl = Rl/cl and fi and a 2 are 
parameters associated with the Roche potential (Webbink 
1977): 

M 2 



Mi + M 2 ' 



a 2 = 



+ 



/<) 



(1 — XLl) 



(11) 

(12) 



where xl,i is the distance of the inner Lagrangian point from 
the centre of the donor in units of a and we have swapped the 
order of the two masses relative to that of Webbink (1977) 
since in our case the secondary star is the donor. 

The key difference between Eqs. |S]and|3]is their sensi- 
tivity to changes in A which has an effect upon when and 
how mass transfer is unstable, but not the rate at which 
mass transfer proceeds when it is stable. Which approxima- 
tion applies depends upon the mass transfer rate, with the 
dividing line given by the mass transfer rate sustainable un- 
der the isothermal model when the inner Lagrangian point is 
located at the photosphere. If the mass transfer rates during 
merging are much higher than this, then the adiabatic rate 
is the more suitable. Following Ritter (1988), the isothermal 
mass transfer rate scales as 

3/2 



M 2 oc RlM^ 1 ( — 



Pph, 



(13) 



where we have left out a weak function of mass ratio. Here 
R 2 , M 2 and T2 are the radius, mass and photospheric tem- 
perature of the donor star while p p h is the photospheric den- 
sity and /1 is the mean molecular mass. For a main-sequence 
star of M 2 = 0.5 M , Ritter (1988) quotes R 2 = 0.52 R Q , 
T 2 = 3520 K, p ph = 1.6 x 10" 5 gcm -3 , p = 1.33 and 
-M 2 = 0.9 x 10 -8 Mq yr" 1 . We scale from this to a white 
dwarf of M 2 = 0.5M Q for which R 2 = 0.013 Rq. Taking 
Koester's (1980) T cS = 20,000 K, logp = 8 model of a DB 
(helium dominated) atmosphere, we find that at the pho- 
tosphere (i = 3.5 and p p h = 2.8 x 10 _5 gcm -3 . From these 
figures we obtain -M 2 = 0.8 x 10~ 12 Mq yr -1 . This rate 
is very much lower than typical equilibrium accretion rates 
at contact which are typically of order 10 -5 Mq yr, except 
in the case of implausibly low component masses. Therefore 
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the adiabatic approximation is the more appropriate. Hence 
for the remainder of this paper, unless stated otherwise, we 
consider the adiabatic model of mass transfer. 

The mass transfer rate contributes towards the evolu- 
tion of the over-fill factor which changes at a rate given by 



It 



where 

C2 = 

Cr L = 



(R.2(2-RL(r 



M2 

M 2 



Rl 



R2 



(C2 



M 2 _a 
Mi a 



•A 



d log i?2 
d log M 2 ' 
d\og(R L /a) 
d log M 2 



(14) 

^W 2 + a ^ 



(16) 
(17) 



The factor Qr L = 1/3 for Paczyriski's (1971) small q approxi- 
mation for the Roche lobe radius while for Eggleton's (1983) 
approximation it is given by 



Cr 



(1 + q) 21n(l + q 1 ' 3 ) - q^/jl + q 1 ' 3 ) 
3 0.6 9 2 / 3 +ln(l + g 1 /3) 



1 



(18) 



for conservative mass transfer, M = 0. White dwarfs typ- 
ically become larger as their mass decreases, and so £2 is 
negative, with typical values of —0.6 < C2 < —0.3. 

Neglecting the second term of Eg. 1151 since |A| <C R 2 , 
and substituting for a /a from Eq.|H] gives 



1 dA 
2R 2 ~~dl 



Jgr _ kMiM 

1 + C2_C^_ 



LO + 



(19) 



■-y/(l + q)r h 



M2 
M 2 



This equation, which is essentially identical to Eq. 20 of 
Priedhorsky & Verbunt (1988), shows that the over-fill grows 
due to gravitational wave losses and shrinks from dissipative 
coupling, if the accretor is spinning faster than the binary. It 
is the last term that can lead to instability, for if the brack- 
eted part is negative then one has a situation where mass 
transfer (M 2 < 0) causes A to grow, which increases the 
mass transfer rate, until the whole process runs away. The 
destabil isation ca used by direct impact is contained in the 
term v/(l"+ q)rh representing the loss of angular momen- 
tum from the orbit. In circumstances that we will elucidate, 
this instability can be counteracted by transferring angular 
momentum back to the orbit via the term in ui. 

In the case of disc accretion, it is usual to assume that 
all the angular momentum contained in the mass gained 
from the donor is fed back to the orb it, which means that 
the terms involving u> and x/(Y+ q)rn do not appear. In this 
case, mass transfer is stable provided that 



q < 1 + 



C2-Cr 



(20) 



because as the mass transfer rate increases and M 2 gets more 
negative, the last term in Eq. 1191 balances the first so that 
A = and M 2 = 0. For a white dwarf donor with £2 ~ — 1/3 
and for C lTL w 1/3, this translates to the widely-used "q < 
2/3" for stability. Relative to this limit, the v/(T~+ q)rh term 
causes instability. However, even with a disc, the material at 
the inner edge of the disc still has some angular momentum. 
In the very tight binary stars we consider here, this turns out 



to be very significant, as we will see in section l2~4l Therefore, 
whether or not accretion is through a disc, we need to include 
the u> term, and therefore to understand how the accretor's 
spin rate evolves with time. 

2.3 The Evolution of the Accretor's Spin Rate 

The accretor will spin up as a result of the addition of mass, 
and spin down because of the dissipative torque, if the ac- 
cretor spins faster than the binary. Conserving angular mo- 
mentum one obtains 



dQ, s 
~~dT 



where 



d0. o dtu 
~~dT + ~dt 



VGMifl fe 
kR\ 



Ma 
Mi 



A _ 1 2 dlog_Rj_ dlogfc 



'dlogAfi d log Mi' 



(21) 



(22) 



a factor which arises from the change in kM^R 1 with the 
addition of mass. 

We mentioned earlier that the moment of inertia factor 
k is approximately 0.2, but more accurately it is a func- 
tion of mass, and decreases at high masses. To allow for 
this we computed values of k as a function of mass from 
Chandrasekhar's equation of state for zero-temperature de- 
generate matter (Chandrasekhar 1967). We fitted these with 
the function 



k = 0.1939(1.44885 - Mi) 



(23) 



which fits to better than 2 per cent over the allowable 
mass range. It should be noted however, that although 
this function is superior to a constant, we are not entirely 
self-consistent since we employ Eggleton's zero tempera- 
ture mass-radius relation, quoted by Verbunt & Rappaport 
(1988): 



R_ 

b' 



= 0.0114 



VMch/ 



-2/3 



M \ 2 / 3 



(m ■,, 



1/2 



(24) 



1 + 3.5 



M 

M v 



-2/3 



+ 



M 
M v 



-2/3 



where M C h = 1.44 M Q and M p = 0.00057 M . The advan- 
tage of this relation is that it provides us with a single re- 
lation that matches Nauenberg's (1972) high-mass relation 
but also allows for the change to a constant density config- 
uration at low masses (Zapolsky & Salpeter 1969). Eggle- 
ton's relation applies for the whole range < M < Mch 
which means that we can use a single relation for both 
white dwarfs, without discontinuities. Our relation for k does 
not account for the effects included by Zapolsky & Salpeter 
(1969), and therefore it is probably an underestimate of k at 
low masses (since k = 0.4 for constant density spheres). The 
denominator in Eq. l24l makes at least a 10 per cent difference 
to the radius for M < 0.07 Mq. However, the lowest mass 
stars that can evolve off the main-sequence within a Hubble 
time have helium core masses of ~ 0.1 Mq, and so it is hard 
to see how white dwarfs of lower mass can be formed. There- 
fore we don't think that the extent of our approximation for 
k will be a problem. 
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Egs 1191 and I21l describe the evolution of the overfill and 
the spin rate of the accretor, but need modification when ac- 
cretion switches from direct impact to disc mode and when 
the donor reaches its break-up rate. We have already men- 
tioned that it is usual to assume that all the angular momen- 
tum is fed back into the orbit when a disc forms (Priedhorsky 
& Verbunt 1988), but because in our case the accretor is rel- 
atively large compared to the binary, in the next section we 
show how to account for the significant angular momentum 
lost at the inner edge of the disc. 



2.4 Disc formation 

As mass transfer proceeds between two white dwarfs, the 
mass ratio become smaller and the orbital separation in- 
creases. Both of these changes mean that there will come a 
point when the stream no longer hits the accretor directly 
and an accretion disc will form. This occurs when the mini- 
mum radius reached by the stream, -R mm , exceeds the radius 
of the primary star, R%. To calculate when this occurs we 
use Eq. 6 of Nelemans et al. (2001). On formation of a disc, 
most of the angular momentum of the accreted material will 
be transferred back to the orbit by tidal forces at the outer 
edge of the disc (Priedhorsky & Verbunt 1988). However, 
material at the inner edge of the disc still has angular mo- 
mentum (e.g. disc accretion acts to spin up the accretor). 
This can be accommodated within Eas ll9l and l21l simmv by 
replacing Rh by R\ and r/, = Rh/a by n = Ri/a whenever 
Ri < Rmin- This is because Rh was defined as the radius of 
a keplerian orbit having the same angular momentum as the 
stream, whereas the inner disc has the angular momentum 
of a keplerian orbit of radius R\. 

Since it is always the case that Rh > -Rmin, when a disc 
starts to form we must have Rh > Ri, so the switch cor- 
responds to a decrease in the angular momentum lost from 
the orbit as a consequence of mass transfer, which therefore 
causes the mass transfer rate to drop. It is important to re- 
alise however that Ri is not that much smaller than Rh and 
so the destabilisation that affects direct impact accretion af- 
fects disc accretion to a very large extent as well. Moreover, 
the square root dependence upon r\ suggests that this effect 
may not be negligible in wider binary systems, a point we 
return to in section |4~S1 

The change from Rh to R\ is the only one we make 
once a disc has formed; we assume that the spin-orbit cou- 
pling expressed through the synchronisation timescale re- 
mains unchanged. 



2.5 Maximum spin rate of the accretor 

For low synchronisation torques, the accretor may be spun 
up to its break-up rate at which point spin-orbit cou- 
pling will presumably strengthen markedly through bar- 
mode type instabilities or the shedding of mass. We model 
this as follows: during evolution if Q s — flk, where 



GMi 
Rf 



(25) 



is the keplerian orbital angular frequency at the surface of 
the accretor, and f2 s > Qk, where 



and 



Ci 



d log Ri 
d log Mi' 



(26) 



(27) 



then we force Cl s = fij. by alteration of the synchronisation 
timescale to a new value t' s given by 



VGWRh \ Mi 
kR\ J Mi 



fit 



(28) 



This revised value is then used in place of uj/t 3 in Eq. ll9l for 
the evolution of the over-fill factor to give 

J_dA _ j GR f | C2~Cr i 

2R 2 dt J oA 



(29) 



M 2 , kMrRfClk 



J orb 



where ri = Ri/a. This modification clamps the spin rate 
of the accretor to the keplerian spin rate at its surface and 
increases the rate at which angular momentum is injected 
back into the orbit. 

Eqs ll9l and l21l together with the above modifications for 
disc formation and the maximum spin of the accretor can be 
integrated to find how the mass transfer rate and the spin 
frequency of the accretor evolve with time. Before looking 
at examples of this, we first look at quasi-static solutions 
of the equations, with our main objective being to uncover 
the conditions under which dissipative coupling can stabilise 
direct impact accretion. 



3 ANALYTIC RESULTS 

The strong dependence of M 2 upon the overfill factor means 
that the timescale of variations in mass transfer is much 
shorter than the mass-transfer (=GR) timescale. Therefore, 
quasi-steady-state solutions, in which the evolution of bi- 
nary parameters is neglected, provide a very useful way to 
understand what happens during mass transfer. To do so we 
assume that the coefficients multiplying uj and M 2 in Eqs ll9l 
and !21l do not change with time and also that Cl — Qk — 0. 
In section [I] we will describe the results of numerical inte- 
grations where we do not make these assumptions. 

The details of the analysis of quasi-static solutions are 
contained in the appendix. We now describe the qualitative 
outcomes of this work, starting with the circumstances un- 
der which dynamical instability is or is not guaranteed to 
occur, and then moving onto the main focus of this paper, 
the stabilising effect of dissipative spin-orbit coupling. 



3.1 Dynamically stable and unstable solutions 

Depending upon the binary parameters, there are three pos- 
sible outcomes: (i) guaranteed dynamical instability, (ii) 
guaranteed stability, and (iii) the intermediate case of ei- 
ther stability or instability, depending upon the degree of 
spin-orbit coupling. It is the third case that is of most in- 
terest in this paper; we leave this to the next sub-section. 
Here we summarise the first two cases. 

Mass transfer is dynamically unstable and will lead to 
a merger, regardless of synchronisation torques, if 
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Figure 1. The upper dashed line shows the dynamical stability 
limit (Eq. 1301 . while the lower dashed line shows Nelemans et 
al.'s stricter criterion (Eg. 1311 . accounting for the switch between 
direct impact and disc accretion at M\ « IMq. The solid line 
shows the transition between disc and direct impact accretion. 
The three dotted lines show how Nelemans et al.'s strict stability 
limit is relaxed when dissipative torques feed angular momentum 
from the accretor back to the orbit (Eq. 1321 , once again account- 
ing for both the direct impact and disc accretion cases. The three 
lines are labelled by the synchronisation timescale in years, the 
line for rg = 1000 yr coincides with the lower dashed line for 
Mi > 1.2 Mq giving the dash-dot line. 
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Figure 2. The evolution of mass transfer rate is shown (solid 
line) for a marginally unstable case of Mi = 0.5 and M2 = 
0.21 Mq and a synchronisation timescale of 30 yr. The model was 
started close to equilibrium and evolved with the masses and or- 
bital separation held fixed. The short vertical lines mark the oscil- 
lation period predicted from a linear stability analysis while the 
dashed curves represent the predicted amplitude. 



typical parameters at birth (section We therefore devote 
a separate section to it, with details once again left to the 
appendix. 



q > 1 + 



C2-C, 



-krj(l + q)\, 



(30) 



(Ea. IA16|l . This is nothing more than the usual "q > 2/3" 
condition for dynamical instability, very slightly modified by 
the last term which accounts for the moment of inertia of 
the accretor. No amount of spin-orbit coupling can stabilise 
the mass transfer in this case. This is the case studied by 
Saio & Jeffery (2002) amongst others; it is also probably the 
most common outcome of a merger as observed double white 
dwarfs have fairly equal masses (Maxted, Marsh & Moran 
2002). 

In the case of direct impact accretion, stable mass trans- 
fer is guaranteed, regardless of spin-orbit coupling, if 



q < 1 + 



C2 - Cr, 



(31) 



(Eq.EIzJ. The same equation with ru replaced by n gives 
the limit in the case of disc-fed accretion. Eq. 1311 is exactly 
that derived by Nelemans et al. (2001) for the case of no feed- 
back of the angular momentum of the accreted material to 
the orbit. The third term on the right-hand side of this equa- 
tion is significantly larger in magnitude than its counterpart 
in Eo. l30l as can be seen in Fig.Qwhich shows parameter con- 
straints for two stars obeying Eggleton's (zero temperature) 
mass-radius relation, Eo. 1241 In this figure, the upper dashed 
line shows the limit, Ea. 1301 for guaranteed instability, while 
the lower dashed line shows the limit, Eq. 1311 below which 
mass transfer is stable, accounting for whether the accretion 
is def by the disc or by direct impact. In between these two 
lines, is a zone in which mass transfer is stable or unsta- 
ble according to the strength of spin-orbit coupling. This 
zone is very significant, both in the fraction of parameter 
space it occupies, and even more so when one accounts for 



3.2 Stability of Equilibrium 

The boundaries of guaranteed stability and instability are 
defined by the existence (or not) of quasi-static solutions of 
Eqs ll9l and l21l The intermediate region of parameter space 
adds the new possibility of unstable quasi-static solutions. 
A linear stability analysis leads to the following condition 
(Eq. IA341 for these solutions to be stable: 



T.S 



> 



P\ 1 + 



Q - y/(l + q)r h + \q 



M 2e 
M 2 : 



(32) 



where — M 2e is the equilibrium mass transfer rate and f3 — 
6R 2 /A e where A e is the equilibrium value of the overfill 
factor corresponding to M 2e (see appendix for details) . is 
a dimensionless factor, typically of order 10 3 to 10 4 , which 
measures the sensitivity of the mass transfer rate to changes 
in the overfill factor. This is where the dependence of the 
mass transfer rate upon A matters: the more sensitive it is, 
the larger /3 is, and, since both the term in brackets multiply- 
ing f3 and M 2e are negative, the stronger the synchronisation 
torque has to be to ensure stability. 

The condition of Eq. 01 is a generalisation of Nelemans 
et al.'s (2001) strict condition for stability (Eo. l3H and is the 
key result of this paper. Once more, this condition applies 
for direct impact accretion, while for disc-fed accretion the 
rh must be replaced by t\ = Ri/a. This condition quanti- 
fies one's expectation that spin-orbit coupling will stabilise 
mass transfer, and essentially says that the synchronisation 
timescale must be less than the timescale upon which the 
mass transfer rate can vary significantly. This is what one 
expects: the spin of the white dwarf must be able to respond 
to variations of the mass transfer rate to ensure stability. 
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Figure 3. The equilibrium spin rate of the accretor Q s rela- 
tive to the keplerian angular velocity at its surface is plotted 
(solid curves) for marginal mass transfer stability for three dif- 
ferent synchronisation timescales (in years). The dashed curves 
show the differential spin rate u> = Q s — Q a relative to Q k . 



(solid line) in Fig. However, at very high accretor masses 
(Mi > 1.2 Mq), even the switch to n is insufficient, and 
the lines of stability drop below the solid line marking the 
disc/direct transition. Thus there are even regimes of disc 
accretion which are destabilised by the loss of angular mo- 
mentum from the inner disc. 

Fig.[3]shows the equilibrium angular velocity of the ac- 
cretor relative to the Keplerian angular velocity at its surface 
in the case of systems just at the stability limit, the fastest 
case. This figure shows that for synchronisation timescales 
of interest for evolution, the accretor does not approach 
break-up. This figure may appear counter-intuitive in that 
weaker synchronisation causes slower rotation in some cases. 
This results from the higher donor masses made possible by 
stronger synchronisation which lead to much smaller orbits 
and orbital periods: for much of Fig. [3] the accretor is al- 
most synchronous with the orbit as can can be seen from 
the small values of the differential spin rate, u = O s — fl 
(dashed lines) , and really we are just seeing that the accretor 
can be close to filling its Roche lobe when it is of low mass, 
and therefore by definition it rotates at a rate of the same 
order of magnitude as the break-up rate. Discontinuities at 
high masses are caused by the transition to disc accretion. 



An example of marginally unstable mass transfer is 
shown in Fig. [5] in which we also compare a numerical inte- 
gration (starting from a slight perturbation of the equilib- 
rium mass transfer rate) with the predictions of the linear 
stability analysis. The oscillations result as first the white 
dwarf is spun up, leading to injection of angular momentum 
back into the orbit which reduces the transfer rate caus- 
ing the white dwarf to spin down, and so on. When the 
synchronisation torque is too weak, the white dwarf does 
not respond fast enough to alterations of the accretion rate 
to damp out perturbations, and their amplitude grows. In 
this particular case the amplitude saturates, but this is not 
of great significance since it is only for rather finely-tuned 
cases that one does not have either stability or such a violent 
instability that merger is inevitable. Moreover, long-term os- 
cillations will not occur in practice because of the evolution 
of the component masses and orbital separation which is not 
included in Fig. [2] 

To apply the stability criterion, Eq. 1321 one must first 
calculate Mie., which also gives A e (through Eq. |UJ; this 
calculation is detailed in the appendix. Lines of stability 
for various synchronisation timescales are plotted in Fig. 
These show how the action of dissipative torques expands 
the region of stable mass transfer in the case of direct im- 
pact. The lower dashed line, marking the onset of instability 
in the absence of any synchronising torque, is raised to be- 
come one of the dotted lines (labelled by the synchronisation 
timescale at the start of mass transfer) . Clearly, the synchro- 
nisation timescale must be short to have much effect, except 
for very low mass systems. Nelemans et al.'s (2001) criterion 
(Eq. l31H is the limiting case for rs — + oo, while the standard 
dynamical stability limit (Eq. 1301 is the limit as rs — > 0. 

For large accretor masses, direct impact can be avoided 
for a wide range of donor masses as the accretor becomes 
very small. The switch from rn to r\ stabilises the mass 
transfer, and for a while stability becomes a case of whether 
the accretion occurs through direct impact or not, with 
the dotted stability lines following the disc accretion limit 



3.3 Super-Eddington accretion 

Eas. lA19l and lA2l"l of the appendix give the equilibrium mass 
transfer rates for dynamically stable systems. Each of these 
is of the form 



-M 2e 
M 2 



/(«)- 



-Jgb. 

"orb 



(33) 



where f(q) — > oo as q increases towards the instability limit. 
Given this, and that for massive systems the GR timescale 
can be short, there are almost inevitably ranges of parameter 
space which, although stable, lead to super-Eddington ac- 
cretion. Han & Webbink (1999) discuss this case extensively, 
focussing upon the fraction of mass that can be ejected. Al- 
though ejection of mass is needed for a system to survive 
super-Eddington accretion, Han & Webbink (1999) argue 
that it is likely that the ejected mass will form a common- 
envelope around the binary system which will lead to the 
merger of the two white dwarfs. Thus we assume, as do 
Nelemans et al. (2001), that the ultimate consequence of 
super-Eddington accretion is merging and the loss of the 
system as a potential AM CVn progenitor. We are in effect 
assuming that dynamical instability and super-Eddington 
accretion lead to the same end result: merging. 

In order to calculate the Eddington accretion rate, we 
follow Han & Webbink (1999) and use the difference in the 
Roche potential at the inner Lagrangian point and the ac- 
cretor, 4>li — 4>a, to give the energy released per unit mass. 
The Eddington accretion rate MEdd is then given by 



A/Edd = 



<JT (4>L1 ~ <f>a, 



(34) 



where gt is the Thomson cross-section of the electron and 
m p the mass of a proton. This expression comes from tak- 
ing two proton masses per free electron, as appropriate for 
fully-ionised helium or carbon. This value is increased over 
the usual solar-composition limit because of the higher mass 
per free electron, but also because of the relatively deep po- 
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tential at the inner Lagrangian point in a double white dwarf 
system. We consider a refinement of this expression in sec- 
tion 



4 NUMERICAL INTEGRATIONS 

The analytical results above take no account of how the sys- 
tem reaches equilibrium or of the evolution of system param- 
eters that occurs during this process. For instance, they take 
no account of the expected lengthening of the synchronisa- 
tion timescale as the binary separation increases, which may 
de-stabilise the mass transfer. Similarly, the analytic results 
do not include the possibility of the accretor reaching its 
break-up spin rate, which increases coupling, and may sta- 
bilise mass transfer. Finally, even though equilibrium mass 
transfer rates can be calculated, it is the maximum mass 
transfer rate that is of more interest from the point of view 
of surviving contact. The significance of these possibilities 
can only be answered through numerical integration. 



4.1 Method 

We carried out fifth-order Runge-Kutta integrations, adapt- 
ing the time-step as the integrations proceeded. The integra- 
tions were started just before contact. To simulate the effect 
of a long interval prior to contact, we started the accretor 
with a differential spin rate of 



dQ 
'~df' 



(35) 



This ensures, through Eq. 1211 that duj/dt ~ in the case 
of strong coupling, as one expects, with the primary star 
lagging slightly behind the increasing orbital frequency. For 
large r s , this can give Q s < 0, in which case we fixed its 
initial value at zero. 

One of the purposes of the numerical integrations was 
to see if systems could be stable and yet exceed the Edding- 
ton luminosity. As described above, Han & Webbink (1999) 
have already detailed the computation of the Eddington rate 
in double white dwarf binary stars. However, they implicitly 
assume that the accretor co-rotates with the binary orbit. 
Since we are explicitly allowing this not to be the case, we 
need to correct Han & Webbink's expression <f>Li — 4>a- In- 
stead we use the following formula 



i>Ll - 4>a — + - (V, - V„) 2 . 



2 -. ■ 2 v- , (36) 

where v< is the impact velocity and v u is the velocity of 
the accretor at the point of impact, both measured in the 
rotating frame. This expression comes from first subtracting 
the contribution of the impact assuming that the accretor 
is co-rotating, and then adding back in the true amount ac- 
counting for uj 7^ 0. If uj = 0, then v„ = and Han & 
Webbink's expression is returned as expected. The modified 
formula works equally for direct impact and disc accretion. 
We implemented it by pre-computing a set of impact veloci- 
ties and locations which we interpolated during the integra- 
tions. In a similar manner we also calculated the luminosity 
of any stream/disc impact in the case of disc accretion. We 
found that once the accretor was spinning at its break-up 
limit, the accretion luminosity was roughly cut in half rel- 
ative to Han & Webbink's (1999) value, which essentially 



is the result of the virial theorem. However, we also found 
that immediately after mass transfer starts, with the accre- 
tor rotating slowly, Han & Webbink's (1999) value is a good 
approximation, and it is only this early phase that matters 
as far as super-Eddington accretion is concerned. 

4.2 Scaling of the synchronisation timescale 

In the analytical calculations, we assumed that the synchro- 
nisation timescale is fixed, whereas one expects it to in- 
crease as the separation of the binary increases. For instance, 
Campbell (1983) calculates the synchronisation due to dissi- 
pation of induced electrical currents within the donor when 
the accretor is magnetic. The synchronisation timescale in 
this case varies with the degree of asynchronism uj, but in 
the small u) limit Campbell (1983) finds 

(37) 



rs 



2 x W 6 M 1 R- 4 R- 5 a 6 B- 2 yv, 



where B is the surface field of the accretor in Gauss and the 
other quantities are in solar units. Putting Mi = 0.6 Mq, 
Ri = 0.012 R Q , R 2 = 0.02 R Q , a = 0.06 R and B = 10 7 G 
gives ts ~ 10 yr. However, this ignores the high conductivity 
of degenerate gas which reduces the dissipation and length- 
ens the timescale considerably for a fixed field strength. For 
example, Webbink & Iben (1987) estimate that the field 
strength would need to be of order 10 10 G to reduce the syn- 
chronisation timescale to the order of decades. One point 
worth making about Campbell's expression is that the a 6 
dependence is largely compensated by the R^ 5 terms since 
J?2 will expand with a, apart from the relatively small ef- 
fect of the changing mass ratio. If this still applies to white 
dwarfs, then it means that if magnetic torques stabilise the 
mass transfer, the systems will remain almost synchronised 
at long orbital periods too. These systems will appear as 
"polars" with magnetically controlled accretion. 

Tides provide another synchronisation mechanism. In 
the standard formalism for tidal synchronisation, tidal de- 
formation in an asynchronous system is damped by some 
form of viscosity or radiative damping (Alexander 1973; 
Zahn 1977; Campbell 1984; Eggleton, Kiseleva & Hut 1998). 
Radiative damping is more effective than the viscosity of de- 
generate matter. Campbell (1984) derives the expression 

'Mi\ 2 ( a \ 6 /M 1 /M N ' 7 



"«» 7 (H) (i) 



yr, 



(38) 



(modified to reflect synchronisation of the accretor rather 
than the donor). Others differ in detail, but retain the scal- 
ing with mass ratio and orbital separation. In contrast to 
the magnetic case, the size of the donor does not enter this 
expression and thus the tidal torque drops off rapidly with 
increasing separation. Unfortunately, as we discuss later, the 
overall magnitude of the timescale seems extremely uncer- 
tain. We therefore take it as a free parameter defined by 
the timescale at the moment of first contact but retain the 
scaling with mass ratio, radius of the accretor and orbital 
separation of the above expression, i.e. we assume that 



rs oc 



(Mi) 2 (— 



We now look at the results of long-term computations of the 
evolution for a wide range of initial masses of each compo- 
nent. 
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Figure 4. The mass transfer rate and accretor's spin rate (relative to the break-up rate) shortly after the start of mass transfer between 
white dwarfs of masses 0.5 and 0.11 Mq with very weak spin-orbit coupling (rg = 10 15 yr). In the left-hand panels the masses were held 
fixed for comparison with analytic predictions (dashed lines). The right-hand panels show that evolution is strong enough in this case to 
avoid the break-up limit, and also to induce a transition to disc-fed accretion after 80,000 yr. 



4.3 No evolution versus evolution of parameters 

Before presenting the full results of the numerical integra- 
tions, we pause to compare a numerical calculation with and 
without evolution (Fig.^J as a test of both the analytic pre- 
dictions (details of which can be found in the appendix) and 
the code. 

This figure shows the switch-on of mass transfer be- 
tween two white dwarfs, for a stable case satisfying Eq. lA17l 
with negligible spin-orbit coupling. Initially the mass trans- 
fer rate rises sharply, but once the mass transfer rate is high 
enough to counter-balance the effect of GR, there follow a 
few thousand years of steady mass transfer at a relatively 
high rate, during which the accretor spins up. The predicted 
rate from the analysis of the appendix is marked as a dashed 
line in the upper- left panel of Fig. 0] and matches the mass 
transfer used by Nelemans et al. (2001) for such cases. Once 
the accretor reaches the break-up rate, then our assumed 
stronger coupling sets in, injecting angular momentum back 
into the orbit and causing a steep drop in mass transfer rate. 

In this particular case (Mi = 0.5 M Q , M 2 = 0.11 M Q ), 
if the parameters are allowed to evolve (as described in the 
next section and plotted in the right-hand panels of Fig.^J, 
then the break-up spin is not reached within the period of 
time shown, but the system makes a transition from direct 
impact to disc-fed accretion which causes a similar, although 
less dramatic, drop in accretion rate after 80,000 yr. 



4.4 Long-term integrations 

With the above scaling of the synchronisation torque, we 
computed evolution over a grid in Mi , M2 parameter space 
to determine the long-term stability of systems. Each model 
was followed for 10 9 yr after contact, or until it became un- 
stable (defined as IM2I > 0.01 M© yr -1 ), or until the ac- 
cretor exceeded the Chandrasekhar limit (defined here as 
1.438 Mq to avoid problems when computing the radius of 
the accretor using Eggleton's formula). We computed the 
grid for a number of different synchronisation timescales at 
contact. The results for ts = 10 15 yr and 10 yr are shown in 
Fig. Qj] while Fig. |S| shows the case of very strong coupling 
with ts = 0.1 yr. In all cases the points plotted represent 
the masses of the stars immediately before the start of mass 
transfer. In these figures we plot the analytic stability lim- 
its as shown in Fig. Furthermore, for the ts = 0.1 and 
ts = 10 15 yr cases we use Eqs. IA19I and IA21I respectively 
to compute the donor mass above which the accretion rate 
will be super-Eddington (an iterative calculation, there is 
no simple functional form for this limit); these limits are 
plotted as dot-dot-dot-dashed lines. 

The main result of our integrations is that neither spin- 
ning up to break-up nor the weakening of the synchroni- 
sation torque as the system evolves have much effect, i.e. 
the analytic stability limit Ea. 1321 provides a fair estimate of 
whether a system will survive mass transfer or not. This is 
because instability either sets in on a very short timescale, 
or not at all. Even those systems which are numerically sta- 
ble while analytically unstable (i.e. they violate the ana- 
lytic stability limit) have little overall effect upon survival 
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Figure 5. The left-panel shows the results of evolution for 10 yr following contact for weak coupling with Tg = 10 yr. The right-hand 
panel shows the case of rg = 10 yr. In each case, models were computed on a regular grid covering the whole parameter range shown, 
but only stable models are plotted, i.e. all the models in the empty upper-left regions are dynamically unstable. Open triangles indicate 
super-Eddington accretion rate systems; open stars indicate super-Eddington accretion rate systems with a total mass in excess of the 
Chandrasekhar limit; filled stars mark systems in which the accretor reaches the Chandrasekhar limit within 10 9 yr. We identify the filled 
circles (stable, sub-Eddington, sub-Chandrasekhar) as AM CVn progenitors. Analytic limits are plotted as in Fig. Qwith the addition 
in the left-hand panel of a dot-dot-dot-dashed line dividing sub- from super-Eddington systems, calculated in the weak coupling limit 
(Eg, IA2TI . Note that this line largely overlaps the dashed line marking instability except for around M\ £S 0.85 Mg and M\ > IMq. 
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Figure 6. As for Figl5l but for very strong coupling (rs = 0.1 yr). 
The dot-dot-dot-dashed line which divides sub- from super- 
Eddington systems were computed in the strong coupling limit 
(Fa, IA19t . The dotted line marking the instability limit for very 
strong coupling almost coincides with the upper dashed line mark- 
ing absolute instability. 

rates because they nearly always exceed the Eddington limit 
(Fig. 0. The small effect that the weakening of the syn- 
chronising torques has can be put down to the concurrent 
lengthening of the GR timescale. We illustrate this in Fig.[7| 
which shows the evolution over 10 10 yr of two weakly coupled 
systems with the same mass for the accretor Mi = 0.5 Mq, 
but with donor masses of 0.125 and 0.126 Mo, the first of 
which is stable in the long term while the second is not. 
This figure shows why our analytic criteria are useful: the 
important action happens so fast that complicating factors 
such as breakup have little effect. 

Three major effects distinguish our numerical integra- 
tions from the analytical work. First, the Eddington lumi- 
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Figure 7. The long-term evolution of two weakly coupled (rg = 
10 15 yr) systems are shown, with one just below the instability 
limit (Aii = 0.5M Q , M 2 = 0.150 M Q , solid line) and the other 
just above (Mi = O.5M , M 2 = 0.151 Mq, dashed). 



nosity is exceeded by some of the dynamically stable sys- 
tems. We assume that these will merge. Super-Eddington 
systems are marked in Figs |5] and |5] by open triangles. The 
second effect, which promotes stability, is the evolution of 
the stellar masses towards the region of stability. This can in 
some cases save a system by moving it from an analytically 
unstable position to a stable one before the instability has 
had time to set in. However, as remarked above, many of 
these "saved" systems suffer super-Eddington accretion, so 
there is little overall increase in survival rates. Third, and 
most obvious, is that some systems are able to exceed the 
Chandrasekhar limit. While it is not clear a priori how long 
they will take to do this, it turns out that the mass transfer 
rate is initially so high, that the total mass of the binary need 
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only be a very small amount over the Chandrasekhar limit 
for this to happen within 10 9 yr, relatively short compared to 
the lifetimes of AM CVn systems. Therefore these systems, 
although of considerable interest as potential Type la su- 
pernova progenitors, will be too short-lived and intrinsically 
rare (because of the large mass of the primary star required) 
to contribute much to the total population of AM CVn sys- 
tems (see Fig. 1101 for example). 

Super- Eddington accretion becomes most significant for 
strong coupling, as seen in the left-hand panel of Fig. |S| In 
this panel, the analytic limit which divides sub- from super- 
Eddington accretion (dot-dot-dot-dashed line, calculated 
from Eq. IA19I agrees well with the numerical integrations. 
However, the analytic limit is obtained in the appendix in 
the limit of rigid coupling between the accretor and the or- 
bit, so no more parameter space can be gained for the pro- 
duction of AM CVn systems by reducing ts any further. 
The analytic limit for the case of weak coupling (Eq. IA2H 
matches the ts = 10 15 yr case well (left panel of Fig. |SJ. 
However, we were not able to obtain accurate analytic pre- 
dictions for the intermediate ts = 10 yr case (right-hand 
panel of Fig. |KJ because in this case it is not correct to as- 
sume either that the accretor is locked as for Eq. IA19I or 
that it is freely rotating as for Eq. IA21I and instead one 
must compute the spin evolution fully. 

We have ignored several effects that may be of signifi- 
cance. First, our models make no allowance for tidal heating 
during the approach to contact (Webbink & Iben 1987; Rieu- 
tord & Bonazzola 1987; Iben, Tutukov & Fedorova 1998). 
This might increase the range of validity of the isothermal 
mass transfer mode, although it is hard to see it having a sig- 
nificant effect given that the transition transfer rate for the 
isothermal versus adiabatic modes is so low. If Campbell's 
(1984) dependence of tidal synchronisation upon luminosity 
applies, there could also be a feedback of tidal heating into 
a stronger synchronisation torque. This could be modelled, 
but with the details of synchronisation so unclear, we prefer 
just to note that such effects might occur. While on this is- 
sue, it is worth emphasizing that we use a zero-temperature 
mass-radius relation for both white dwarfs, whereas we can 
expect larger radii for a given mass from finite entropy ef- 
fects, depending upon the time taken to reach contact (De- 
loye & Bildsten 2003). It is also worth pointing out that 
if tidal heating is not significant, then the donor may not 
even be synchronised, which would necessitate changing the 
whole prescription based upon Roche geometry. Another is- 
sue that we have not tackled is the effect of helium ignition 
on the accreting white dwarf (since the donors of interest 
for survival of mergers are all low mass and therefore of he- 
lium rather than carbon-oxygen composition). It is not clear 
whether helium ignition, if it takes place, would increase or 
decrease the probability of avoiding the merger. 

4.5 Destabilisation of Disc Accretion 

It is the finite size of the accretor which leads to the destabil- 
ising removal of angular momentum from the orbit as mass is 
transferred. Direct impact accretion is a secondary, although 
important, consequence of this. To illustrate this point, we 
repeated the weak coupling calculations shown in the left- 
hand panel of Fig. while assuming that all accretion oc- 
curs through a disc, regardless of the whether the minimum 




Mj (M ) 

Figure 8. This figure again shows very weak coupling (ts = 
10 yr, cf left-hand panel of Fig[5J, but with the accretion forced 
to be through a disc. The lower dashed line shows the guaranteed 
stability limit for disc accretion, while the thin solid line shows 
the usual limit when direct impact is included. 

stream radius was smaller than the primary star or not (one 
could envisage for instance initialising these systems in a 
state of disc accretion, however unrealistic this is in prac- 
tice). The results are shown in Fig. [H] Rather surprisingly 
perhaps, we see that the standard "q = 2/3" limit, marked 
by the upper dashed line, is not much better for disc accre- 
tion than it is for direct impact accretion. This is also clear 
from the relatively small change between the thin solid line 
(direct impact included) and the lower dashed line (disc ac- 
cretion only). As we said earlier, this is because the primary 
star's radius Ri is not that much smaller than the circular- 
isation radius Rh in these systems. For disc accretion, the 
zero spin-orbit coupling stability limit of Eq. l31l becomes 



g <l + ^_^-^(l + g )^, (40) 

which is plotted as the lower dashed line in Fig. |H] 

This raises the question of whether sinking angular mo- 
mentum into the accretor also has a significant effect in cata- 
clysmic variable stars (CVs). To evaluate this, we computed 
the analytic stability limit with a revised mass-radius rela- 
tion for the donor. We assumed that R2 = M2 (solar units) 
and that £2 = —1/3 (measuring the adiabatic response), ap- 
proximately correct for low mass donors. The result is shown 
in Fig. |U] As expected, given the larger orbits of CVs, the 
destabilisation is nothing like as severe as it is for double 
white dwarfs, but neither is it entirely negligible. The im- 
pact of this upon the formation of CVs depends upon the 
distribution of pre-CVs in the (M\,M-i) plane, but has some 
potential to favour magnetic over non-magnetic systems and 
to reduce overall formation rates. 



5 SURVIVAL AS A FUNCTION OF 
SYNCHRONISATION 

We are unlikely to catch a system in the short-lived phase 
after the start of mass-transfer, and so the most important 
consequence of the destabilising effect of finite accretor size 
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Figure 9. The regions of stability for cataclysmic variable stars 
with negligible spin-orbit coupling (t<j = 10 yr). Systems below 
the dotted line (Eq. gHJ are stable, whereas under the usual as- 
sumption that angular momentum is entirely fed back, the upper 
dashed line would be the relevant limit. The lines on this figure 
are the CV equivalents of those plotted in Fig. [S] 




Figure 10. The masses of double white dwarfs at birth from the 
binary population models of Nelemans et al. (2001) that survive 
in the case of strong coupling (t$ = 0). The surviving models 
do not reach the limit of guaranteed stability in this case (upper 
dashed curve) because super-Eddington accretion becomes the 
more stringent constraint (as marked by the dot-dashed line). 
Dotted lines are as plotted in Fig. Q 



is its impact upon the survival of double white dwarfs as 
AM CVn systems. To gain an idea of how significant this 
could be, we applied our stability criterion to the models of 
Nelemans et al. (2001) to obtain the predicted birth rate of 
AM CVn stars as a function of the synchronisation timescale 
at the start of mass transfer. Fig. llOl shows a greyscale repre- 
sentation of the stellar masses at birth from the models. The 
greyscale only shows systems that can become AM CVn bi- 
naries. They concentrate towards the line dividing sub- from 
super-Eddington systems because they are outliers from the 
majority of systems which have near-equal mass ratios. This 
depends upon the physics of the common envelope amongst 
other things: Nelemans et al. (2001) developed their model 
partly in response to the observations of equal mass ratios 




Figure 11. The birth rate of AM CVn systems descended from 
double white dwarfs as a function of the synchronisation timescale 
at the start of mass transfer. We show also the case for isothermal 
mass transfer (dotted line) to show that the mass transfer mode 
has a relatively slight overall effect. The plot was derived from 
our analytic approximations for speed, but limited comparisons 
with numerical integrations produced nearly-identical results. 



of double white dwarfs. Others have found more unequal 
distributions (Iben et al. 1997; Han 1998), nevertheless all 
models so far peak with q > 0.4 and thus all predict that 
most systems are unstable and will merge. The high den- 
sity near the stability limit of Nelemans et al.'s distribution 
means that the birth rate of AM CVn systems from the 
double white dwarf path is very sensitive to the strength of 
synchronisation torque. This can be seen in Fig. llll in which 
we plot the birth-rate (from the double white dwarf merger 
route) as a function of the synchronisation timescale 1 . This 
figure shows that as spin-orbit coupling weakens, the birth 
rate of AM CVn stars from double white dwarfs (which we 
assume to be equivalent to survival of mass transfer as a bi- 
nary) drops by about one hundred-fold. A synchronisation 
timescale ts < 1000 yr is enough to increase the birth rate 
from its lowest value by more than a factor of two. While 
sensitive to the mass ratio distribution of the detached dou- 
ble white dwarfs, this calculation illustrates the potential 
significance of dissipative synchronisation torques for the 
population of AM CVn systems. In this figure we also show 
the results for isothermal mass transfer, which shows that 
variations between the two extremes is not that large. 

One could hope that an observational estimate of the 
birth rate of AM CVn stars might decide whether we are in 
the strong or weak coupling part of Fig. 1111 Unfortunately in 
our opinion, uncertainties and selection effects are too large 
to allow this to be done. First, there are only two directly- 
measured distances (CP Com, 68±7pc, Thorstensen 2004; 
AM CVn, 235 pc, C.Dahn 2 , priv. comm.). The AM CVn dis- 
tance is a factor of 3 further than Warner (1995) assumed in 
estimating the space density, which if it applied to all sys- 



1 The extreme values in Fig. llll diffcr from those of models I and 
II in Nelemans et al. (2001) because we use a different mass-radius 
relation in the calculation of the Eddington rate. 

2 On behalf of the USNO CCD parallax team 
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tems would imply an overestimate by a factor of 27, which 
serves to show how uncertain observational estimates are. 
To add to this, Nelemans et al. (2001) estimate that while 
99% of AM CVns have periods P > 2000 sec, they will only 
form ~ 10% of the observed sample, an estimate that itself 
depends upon highly uncertain assumptions about selection. 
Finally, the models that lead to Fig. llll have uncertain nor- 
malisation - a factor of 5 to 10 seems quite possible. There 
is clearly considerable scope for work in this area. 



6 DISCUSSION 

The requirement that the initial synchronisation timescale is 
below 1000 yr in order to raise survival rates of double white 
dwarfs as binary systems is the most important result of 
this paper. There are few calculations of synchronisation of 
white dwarfs in binaries, but those that exist come nowhere 
near this strength of synchronisation. For instance, Camp- 
bell's (1984) estimate for radiative damping (Eq. 1381 ap- 
plied to a pair of white dwarfs with masses of Mi = 0.5 Mq 
and M 2 = 0.2 Mq gives ts ~ 10 12 yr. Even longer times 
(10 15 yr, Webbink & Iben 1987) emerge from calculations 
based upon the viscosity of degenerate matter. Indeed, the 
estimates of synchronisation timescales are so long that even 
the mass donor would not have achieved synchronism prior 
to mass transfer let alone the accretor (Mochkovitch & Livio 
1989). If these estimates are reliable, then the lowest birth 
rate estimates of Nelemans et al. (2001) are to be preferred. 
However, it is widely recognised by the authors of the pa- 
pers on synchronisation of white dwarfs that the problem 
is essentially unsolved: the possibilities of turbulent viscos- 
ity (Horedt 1975) and the excitation of non-radial modes 
(Campbell 1984) can dramatically shorten the synchronisa- 
tion time. For instance, Mochkovitch & Livio (1989) argue 
that turbulent viscosity can give synchronisation timescales 
<C 500 yr, exactly what is needed to raise the birth rates sig- 
nificantly. The uncertainty over the dissipation mechanisms 
within white dwarfs is thus a major unsolved problem for 
AM CVn evolution. 

If the synchronisation torques are not strong, then ei- 
ther the ultra-short period systems (if they are such) are de- 
scended from systems that have always had extreme enough 
mass ratios to be stable, or the double white dwarf route may 
not be the main channel for these systems. The former alter- 
native is possible, despite the absence of a single system of 
extreme mass ratio amongst the observed close double white 
dwarf population (Maxted et al. 2002) , because they are ex- 
pected to be rare, as shown by Fig. 1101 This would mean 
that AM CVn birth rates through the double white dwarf 
route are dependent upon on a so-far-untested part of the 
mass ratio distribution of double white dwarfs. The main 
objection to formation other than from two white dwarfs 
is the relative difficulty of reaching ultra-short orbital peri- 
ods by other routes: it is important from this point of view 
that the photometric periods of order 10 minutes or less are 
confirmed (or not) as orbital periods. 

One way to test our ideas would be to establish that 
there are systems in existence that are stabilised by syn- 
chronisation torques. Although such systems are likely to 
be in a phase of direct impact accretion, it is important to 
realise that direct impact on its own does not equal instabil- 



ity. There is a region of parameter space of stable, direct im- 
pact accretion even in the absence of synchronisation torques 
(the slim lenticular region in Fig.Qcontaining the word "Di- 
rect"). For instance, V407 Vul, if it is a direct impact ac- 
cretor, could be in this region (Marsh & Steeghs 2002); so 
too could RXJ J0806. 3+1527. Thus even if direct impact 
accretion were to be established in a double white binary, it 
would not mean that significant synchronising torques were 
acting. We would need in addition to establish that the sys- 
tem parameters implied instability in the absence of such 
torques. Short period systems are of particular interest in 
this respect since short periods require high donor masses, 
and correspondingly high accretor masses if they are to re- 
main in the lenticular region. This is an excellent reason 
for attempting to find further examples of such systems. X- 
ray surveys would seem the most promising method, since 
both V407 Vul and RX J0806.3+1527 were found in this 
manner, although in the future perhaps space-based gravita- 
tional wave detectors will prove even more useful (Nelemans, 
Yungelson & Portegies Zwart 2004). 

Similar remarks apply to CVs, except for these stars the 
unstable region of Fig. [5] is small enough that it seems un- 
likely that the stellar masses can be measured with enough 
accuracy to be certain that a given system is located within 
it. 

6.1 Measuring synchronisation torques from other 
binary stars 

The strength of synchronisation torques on white dwarfs 
might be testable in other systems. If any short period sys- 
tems with rapidly rotating white dwarfs can be identified, 
it might be possible to place a lower limit upon the syn- 
chronisation torque acting. For instance, rotation velocities 
have been measured for several white dwarfs in CVs (Sion 
1999). The strongest constraints come from high spin-rate, 
low accretion-rate, short-period systems. Perhaps the best 
example is the cataclysmic variable WZ Sge, which has an 
orbital period of 81.6 min and an accretor which shows signs 
of having a spin period of « 28 sec (Patterson et al. 1998). 
WZ Sge is an old, low accretion rate system, accreting via 
a disc. Assuming that its spin period represents the equilib- 
rium between accretion and tidal dissipation torques, and 
since Q s 3> Q , we have 

r s = — , (41) 

\m 2 \VgWK 

where R a is the radius equivalent to the specific angular mo- 
mentum of the accreted material. We expect R a to lie in the 
range Ri < R a < Rh (allowing R a ^ Ri because the accre- 
tor may be magnetised given the presence of a spin signal) . 
We take for the mass transfer rate \M 2 \ = 3x 10 -11 Mq yr~\ 
estimated from Fig. 1 of Kolb & Baraffe (1999) as we are 
interested in the long term average value since the timescale 
turns out to be long. Thus assuming for WZ Sge that 
Mi = 1M Q and M 2 = 0.08 M Q (Steeghs et al. 2001), we 
find 

3 x 10 8 < T S < 1.4 x 10 9 yr. (42) 

Armed with the constraint on ts from WZ Sge, we can 
scale to a double white dwarf with Mi = 1 Mq and M 2 = 
0.26 Mq . The separation of such a system is 9.5 times smaller 
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than WZ Sge, while the donor is 3.25 times more massive. 
Scaling according to ts oc a e M r2 (Mi and Ri being the 
same) , the equivalent range of ts for the double white dwarf 
is 

40 <t S < 180 yr, (43) 

very much the sort of values that can have an impact upon 
AM CVn evolution. 

Unfortunately, there are significant problems with the 
above analysis. If the pulsation in WZ Sge really does mean 
that the accretor is magnetic (and truly represents the white 
dwarf's spin period), then it may be magnetic coupling to 
the outer accretion disc which keeps it rotating slowly. In 
addition there are complications of nova outbursts in CVs 
which probably make them unsuitable for measurement of 
ts under any circumstances (Livio & Pringle 1998). De- 
tached systems avoid these problems, and perhaps might 
make it possible to measure a spin-down rate if a weakly- 
magnetic example could be found (but strong enough to 
give a spin pulse) and it could be shown that magnetic and 
accretion torques were negligible. Interestingly there is ap- 
parently a detached system which contains a fairly rapid ro- 
tator (EC13471-1258, O'Donoghue et al. 2003), but in this 
case it is suggested that this was caused by a recent episode 
of accretion, and so it is not an indication of weak synchroni- 
sation. We are left with no clear answer as to the strength of 
synchronising torques from the known types of white dwarf 
binary stars. Nevertheless, it is clear from the estimate above 
that if tidal synchronisation is significant for the evolution 
of AM CVn systems, then it should also be important in 
fixing the spin rates of white dwarfs in cataclysmic variable 
stars and similar systems, and this could be a possible ex- 
planation for the generally low rotational velocities detected 
so far in cataclysmic variable stars, none of which are near 
the break-up rate (Sion 1999). 

6.2 Magnetic versus tidal synchronisation 

If magnetic rather than tidal torques were dominant, an 
interesting filtering effect could occur in which systems 
with strongly magnetic accretors survive mass transfer while 
non-magnetic systems mostly merge. One would then ex- 
pect magnetic accretors to be over-represented amongst 
the AM CVn systems. This is not obviously the case, e.g. 
CP Com has no detectable circular polarisation (Cropper 
1986) and along with V396 Hya shows the clear signature 
of emission from a disc, while other systems in the class 
show shallow absorption line spectra similar to those of high- 
state non-magnetic CVs rather than the magnetic AM Her 
class. This suggests that either the fraction of merging dou- 
ble white dwarfs that are magnetic is so low that even the 
pruning of 99 per cent of the non-magnetic systems fails to 
reveal them, or that tidal synchronisation is effective in sta- 
bilising both non-magnetic and magnetic systems, or that 
neither magnetic nor tidal synchronisation has much effect 
or, finally, that AM CVn stars do not form through the dou- 
ble white dwarf route. 

6.3 Period changes 

As several investigators have realised, measurements of pe- 
riod changes have some potential to discriminate between 
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Table 1. The ratio of the period derivative in a semi-detached 
state to its value when (just) detached is given for different stellar 
masses. 

models of the ultra-compact binary stars. Once a sys- 
tem has settled onto its long term evolution towards the 
AM CVn phase and spin equilibrium has been reached, it 
is straight-forward to show that the period derivative in a 
semi-detached state Psd relative to its pure GR-driven value 
when detached Pr>, is given by 

Pan = (Cs-O-J/2 m 

Since £ r/ . ~ 1/3 and £2 ~ —1/3, and since the denominator 
of the above expression must be positive for stability, this 
ratio is negative; see Table^for some explicit values for sys- 
tems which are stable even when there are no synchronisa- 
tion torques. To date period changes have been reported for 
V407 Vul (Strohmayer 2002) and RX J0806.3+1527 (Hakala 
et al. 2003; Strohmayer 2003), and in each case the period 
was found to be decreasing. These results favour a detached 
rather than semi-detached configuration for these two sys- 
tems (ES Cet on the other hand is clearly accreting). This 
would weaken the case for V407 Vul and RX J0806. 3+1527 
having passed through the state discussed in this paper, al- 
though period changes in other close binary stars such as 
the cataclysmic variable stars have proved unreliable as a 
means of measuring long-term angular momentum changes 
(Applegate 1992; Baptista et al. 2003). The periods of the 
ultra-compact systems should continue to be monitored to 
see whether similar problems affect these stars, although one 
might hope that at such short periods gravitational wave 
angular momentum losses will dominate over any secondary 
effects. 



7 CONCLUSIONS 

We have studied the onset of mass transfer between two 
white dwarfs, in particular its stability during the phase 
when the accretion stream hits the accretor directly. We 
find that this phase can be stabilised by coupling of the 
spin of the accretor to the binary orbit, through dissipative 
processes such as tidal stressing and magnetic induction. 
However, the coupling needs to be strong and must act on a 
timescale < 1000 yr when mass transfer first starts to have 
much effect upon the survival rates of the systems. Standard 
estimates of the synchronisation timescales in white dwarfs 
are orders of magnitude longer than this, but may wildly un- 
derestimate the synchronisation torques. We have also found 
that during disc accretion, the angular momentum lost to 
the accretor at the inner edge of disc also destabilises mass 
transfer, and is almost as significant as the direct impact 
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case, essentially because the radius of the accretor is com- 
parable to the circularisation radius in these systems. The 
same effect plays a lesser but non-negligible role in non- 
magnetic cataclysmic variable stars. Further discoveries of 
ultra-compact double white dwarfs have the potential to tell 
us whether synchronisation torques are indeed important for 
surviving this exciting phase of evolution. 
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APPENDIX A: QUASI-STATIC SOLUTIONS 

In this appendix we are interested in quasi-static solutions of 
equations 1191 and 1211 "static" on timescales short compared 
to the secular evolution of the system parameters and their 
stability. As outlined in section Um we assume therefore that 
the coefficients involving binary masses and separation are 
constant. 

We define the following dimensionless variables 
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the first being a rescaling of the timescale. With these new 
variables, equations 1191 and 1211 can be written as 



dz 

— = oi — a 2 y — 032; 
at 



dt 



-a^xy — y + a$x, 



(A5) 
(A6) 



where the dimensionless coefficients a± to as are given by 
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(A9) 
(A10) 

(All) 



and where tq = — Jorb/^GR and rj = 7?i/a as before. The 
coefficients ai, a 2 and 05 are positive, while 03 and 04 can 
be either positive or negative. For Eggleton's mass-radius 
relation, Eq.|^] for instance, a 4 < for Mi > 0.40 M Q . 

For the solutions of interest, y = z = 0, and we are left 
with two simultaneous equations 



a 3 x e + a 2 y e = ai, 
a5X e —y e = a 4 x e y e . 

These lead to a quadratic equation with solutions 



—b ± \/b 2 + 4aia3<i4 
20304 



where 



b = as + a 2 a 5 — a\a 4 . 



(A12) 
(A13) 



(A14) 



(A15) 



Solutions of physical interest must have real, positive x e , 
and, for reasons of stability, if there are two such solutions, 
the smaller of the two is favoured. 

This divides parameter space into two parts according 
to whether there are or are not any solutions. The region 
where there are solutions is further sub-divided according 
to whether or not the solutions are stable. 



AI The dynamically unstable case 

We know that ai > 0. Hence, if 04 > and 03 + a 2 as < 
(which since a 2 and as are both positive implies as < 0), 
then b < and \b 2 + 4ctia3a4| < \b\, and so there are no real 
and positive solutions for x e . The condition as + a 2 a$ < 
is equivalent to: 



q > 1 + 



krt(l + q)\, 



(A16) 



which was quoted in the main text as Eq. 1301 

If 04 < 0, as is possible, there is formally a positive 
solution of Eq, IA14l but only at such high spin rates that we 
run into the break-up limit of the accretor. The consequence 
is that Eq. lATol remains a good condition for instability. 



A2 The dynamically stable case 

The reverse of the unstable case is when as > 0, when one 
can show that there is always a solution for real, positive 
x e , for any ts and a 4 . The condition 03 > means that 



q<l 



+ q)r h , 



(A17) 



as derived by Nelemans et al. (2001). As noted in the main 
text, this can be adapted for disc-fed accretion by replacing 
r h by n = Ri/a. 

Systems evolve towards this solution with time since 
their mass ratios continuously decrease. We will prove below 
that the equilibrium is stable in this case. Before doing so, we 
look at the equilibrium mass transfer rate in this stable case. 
The equations derived here are those used in the comparison 
of analytical and numerical results shown in Fig. 2] 

There are two extreme cases as ts — > or ts — » 00. As 
ts —> 0, ai — > and 

x e -> ^ . (A18) 

a 3 + a 2 a 5 

Dividing out the ts in the definitions of x and Oi, then we 
have 



-M 2e 
M 2 



— Jgr/Jo: 



1 +(Ca-Gx)/2 ■ 



krl\(l + q) 



(A19) 



This is, as expected, the usual expression for stable mass 
transfer, albeit slightly modified by the term including kr\ 
as a result of the angular momentum removed from the orbit 
by the accreting white dwarf owing to its increasing mass. 
This term contributes at most ~ 0.01, and can safely be 
ignored. 

When ts — > 00, we run into break-up of the accretor 
once more. It is then clear from Eq. l29l that the equilibrium 
mass transfer rate is given by 



-M 2e 
M 2 



-Jgr/ Jorb 



1+ (Ca - CrJ/2 - q - k\y/(l + q) ri 



(A20) 



The final factor in the denominator is again fairly small, and 
so this mass transfer rate is very similar to that of Eq. IA19I 
Thus we conclude that systems undergoing direct impact ac- 
cretion which satisfy Nelemans et al.'s (2001) strict stability 
condition, Eq. 1311 always have equilibria with mass transfer 
rates close to that of Eq. IA19I 

Nelemans et al. (2001) use a different equation to com- 
pute the Eddington-limited mass transfer rate in the stable 
direct impact case. On their Fig. 1 they plot a dashed line 
obtained from the relation 



-M 2e 
M 2 



1+ (Ca - CrJ/2 - q - ^(l + q)r h 



(A21) 



which gives substantially higher mass transfer rates than 
either of Eos lA19l and lA20l Although this is not the equilib- 
rium rate, Nelemans et al. were in fact correct to use this, 
at least in the weak coupling case, as can be seen by consid- 
ering the situation that applies as ts — > 00 while to remains 
finite. Eq. I19l then becomes 



1 dA 
2R 2 H 



Jgr 

J orb 



+ 1 



\/{l+q)r h 



M2 

M 2 



(A22) 



and the equilibrium solution, A = 0, gives Eq. IA21I Al- 
though this a temporary state that will stop once the white 
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dwarf has spun up to break-up, it is still long-lasting (see 
Fig. [*"J. The spin-up phase will last a time tb given by 



ky/n/r h (l + (Ca - Cr J/2 - « - V(T+ ?)r h ), (A23) 



Tb 
tg 

which can be long enough for the system to switch to disc ac- 
cretion, although once break-up is reached, the mass transfer 
rate is still given by Eq, IA20l These results are the ones illus- 
trated graphically in Fig. [I] As spin-orbit coupling increases 
we expect a peak rate intermediate between that given by 
Ea. lAlOl and that given by Eq. lXSTI 



A3 Stability of Equilibrium 

The case of most interest to us is when the mass ratio lies 
in between the unstable and stable limits of Eqs IA16I and 
IA17I when —a 2 a^ < 03 < 0. In this case synchronisation 
can bring about an equilibrium solution which disappears 
as ts — ► 00. This then sets an upper limit on ts required for 
the existence of an equilibrium solution. However it turns 
out that an even stricter constraint comes from the require- 
ment that the equilibrium solution be stable. This is a new 
feature: so far there either have been or there have not been 
equilibrium solutions, and if they existed we have implic- 
itly assumed that they were stable equilibria. We now come 
across equilibrium solutions that can be stable or unstable. 

Consider the following small perturbations from equi- 
librium: 



V = Ve+y , 
Z = z e +z, 

X — Xe. ~\~ (3Xe% j 



(A24) 
(A25) 
(A26) 



where in the last equation we recognize that x varies with 
z alone since the mass transfer rate depends explicitly upon 
A but not (j, and the dimensionless factor p is defined by 



P = 2R 2 



dlog(-M 2 
dA 



(A27) 



Since the mass transfer rate increases monotonically with 
A, p is positive. For the two mass loss models we consider 
in section 12*751 

Pi = 2^, (A28) 
for the isothermal case of Eq. [HI and 

Pa = 6 J, (A29) 

for the adiabatic case of Eq. |"3 where A e is the equilibrium 
value of the overfill factor; it is the latter expression that we 
employed throughout the paper. 

Applying the perturbation to Eqs IA5I and IA6I while 
keeping terms to first order in the small dashed quantities 
gives 



~f3x e a,3z' — a,2y . 
f3x e (a,5 — a4y e )z' 



(1 + a 4 x e )y 



(A30) 
(A31) 



These are two first-order, coupled equations. Assuming solu- 
tions of the form e pt , one obtains an eigenvalue equation for 
p with values that are solutions of the following quadratic: 



after using Eqs lA12l and lA13l to simplify the third term. The 
solutions are stable provided that Re(p) < for both roots. 
This is true if the two coefficients of the quadratic are both 
positive. The coefficient of the term linear in p is positive if 



— > - [a s f3 + a 4 ] - 1 r j —, 
T S M 2 

or, when written out in full, 



(A33) 



ts 



1 + 



+ q\ 



-Ma. 
M 2 



V(l + g)r h ^ P (A34) 



which was quoted in the main text fEa.!32[l. It turns out that 
the coefficient of the constant term is also positive when this 
condition is satisfied, so no further restriction is necessary. 
Eo, IA34l is a generalisation of Nelemans et al.'s (2001) strict 
condition for stability C Eg. 1311 and is the key result of this 
paper. Eq. IA34I derived from the condition, a^P + 0,4 > 0, 
does not quite match the criterion for stability regardless of 
ts, Eq. IA171 which corresponds to as > 0. This is because 
Eo. IA341 permits arbitrarily high rates of spin, whereas in 
reality the spin rate is limited by break-up. Both criteria 
are correct, but under different conditions: Eo. IA17I applies 
in the case of zero synchronisation, while Eq. IA34I applies 
when there are synchronising torques strong enough to keep 
the spin rate of the accretor below the break-up rate. 



p + (1 + a4X e + Px e as)p + P (ai + a^a^x,. 



0. 



(A32) 



